No separate roles of repulsive and attractive pair forces for the physics of liquids 
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According to standard liquid-state theory repulsive and attractive pair forces play distinct roles 
for the physics of liquids. We here present a continuous series of pair potentials that have virtually 
the same structure and dynamics, although only some of them have attractive forces. This reflects 
the fact that the motion of a given particle is determined by the total force on it, not the individual 
contributions to the force from other particles, which in terms of the pair potential is the quantity 
usually discussed in liquid-state theory. 



A liquid is held together by attractions between its 
molecules. It is very difficult to compress a liquid be- 
cause the molecules strongly resent closely approaching 
each other. These facts have been well known for a long 
time, and today it is textbook knowledge that the repul- 
sive and the attractive forces play distinct roles for the 
physics of liquids. The repulsive forces, which ultimately 
derive from the Fermi statistics of electrons, are harsh 
and short ranged. According to standard theory these 
forces are responsible for the structure and, in particular, 
for reducing considerably the liquid's entropy compared 
to that of an ideal gas at the same density and tem- 
perature. The attractive forces, on the other hand, are 
long ranged and weaker. The attractions, which derive 
from induced dipolar interactions, reduce the pressure 
and energy compared to that of an ideal gas at the same 
density and temperature. We argue below that this pic- 
ture, though very appealing, overemphasizes the individ- 
ual pair forces and does not provide a full understanding, 
because systems may be constructed where the repulsive 
and attractive interactions clearly do not play separate 
roles. 
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FIG. 1: Lennard- Jones pair potentials VLj(r) = ie[(r /a)~ 12 — 
(r/(j) -6 ] predicted to give the same physics at the state point 
(p,T) = (1,1). Visually these potentials have little in com- 
mon; in particular, they have very different contributions from 
attractive forces. The pair potentials were constructed ana- 
lytically using the isomorph theory (detailed in the text after 
Fig. 0. 
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The traditional understanding of the liquid state is 
based on pioneering works by Frenkel, Longuet-Higgens 
and Widom, Barker and Henderson, and Weeks, Chan- 
dler, and Andersen (WCA), and many others fH, 12]. The 
fundamental idea is that the attractions may be regarded 
as a perturbation of a Hamiltonian based on the repulsive 
forces, the physics of which is usually well represented by 
a hard-sphere reference system [3| ■ Perturbation theories 
based on this picture are standard for calculating 

simple liquids' thermodynamics and structure as quan- 
tified by the radial pair distribution function g(r). We 
do not question these very useful perturbation theories, 
but argue below from theory and simulations that the re- 
pulsive and the attractive pair forces do not always play 
clearly distinguishable roles for the structure and dynam- 
ics of simple liquids. 



This point is illustrated in the simplest possible way by 
studying systems of Lennard- Jones (LJ) particles. The 
LJ pair potential is given by ULj(f) = 4e[(r/o-)~ 12 — 
(r/er)~ 6 ], where e and a specify the energy and length 
scales, respectively. It is shown below that the different 
LJ pair potentials of Fig. [T]give rise to virtually the same 
structure and dynamics at the state point (p,T) — (1,1) 
(this defines the unit system used below and is the unit 
system in which the e's and cr's of Fig. [T]are given). The 
potentials of Fig. Q] all have attractive forces, but for 
some of the potentials the attractive forces are insignifi- 
cant. To show that these potentials have the same struc- 
ture and dynamics, NVT computer simulations of sys- 
tems of 1,000 particles were performed using our RUMD 
software, which runs on graphical processing units 
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FIG. 2: (a) The radial pair distribution function for the LJ 
pair potentials of Fig. [T] at the state point (p,T) = (1,1) 
(left panel) and for a series of LJ pair potentials with fixed 
a parameter (right panel) at the same state point, (b) The 
mean-square displacement for the LJ pair potentials of Fig. [1] 
at the state point (p,T) = (1, 1). 



Figure (5Ja) (left panel) shows the radial distribution 
function g(r) for the seven LJ pair potentials of Fig. Q] 
at the state point {p,T) = (1, 1). For comparison, simu- 
lations are shown in the right panel for seven potentials 
with the same e variation, but fixed a = 0.947. Figure 
[2Kb) shows the mean-square displacement for the Fig. [1] 
potentials. By the Henderson uniqueness theorem Q the 
pair potentials of Fig. [1] cannot have exactly the same 
pair distribution functions, but clearly they do have lead 
to very similar structure and dynamics. 

How were the pair potentials of Fig. [T] determined, 
and why do they have the same structure and dynam- 
ics? The starting point is the existence of isomorphs in 
the phase diagram of liquids with strong correlations be- 
tween NVT virial and potential-energy equilibrium fluc- 
tuations @, Q (which we recently argued is a useful def- 
inition of a "simple" liquid Q). Two state points with 
density and temperature (pi,Ti) and (p2,?2) are termed 
isomorphic [7fl if all pairs of physically relevant micro- 



configurations of the two state points, which trivially 
scale into one another [i.e., p± T^ — p\ T? for all 
particles i] , have proportional configurational Boltzmann 
factors. This means that exp[— U(r^\ ...,r^J')/fcsTi] = 



(2,) [ ' i 

C12 exp[— U(r\ , r^ ; )//cBT 2 ]. L J systems are strongly 
correlating and thus have isomorphs to a good approxi- 
mation [8j . The invariance of the canonical probabilities 
of scaled configurations along an isomorph has several 
implications [7J: Excess entropy and (excess) isochoric 
specific heat are both isomorph invariant, the dynamics 
in reduced units are invariant for both Newtonian and 
Brownian equations of motion, reduced-unit static den- 
sity correlation functions are invariant, a jump between 
two isomorphic state points takes the system instanta- 
neously to equilibrium, etc. For Newtonian dynamics 
using reduced units corresponds to measuring length in 
units of p -1 / 3 , time in units of p~ x ^ ^mfk^T where m is 
the particle mass, and energy in units of fcgT. Thus the 
reduced particle coordinates are defined by r* = p 1 / 3 !^. 

Recently it was shown that liquids with good iso- 
morphs have simple thermodynamics in the sense that 
the temperature is a product of a function of excess en- 
tropy per particle s and a function of density [loj |. 



.(2) 



T = f(s)h(p). 



(1) 



The function h(p) inherits the analytical structure of the 
pair potential in the sense that if the latter is given by the 



expression v(r) = J3 T 



then h(p) = J2n °nP n/3 



in which each term corresponds to a term in the pair 
potential. Since h(p) is only defined within an overall 
multiplicative constant, one can write for the LJ pair 
potential 



h(p) = ap A + (1 - a)p 2 



(2) 



The constant a was determined from simulations at the 
reference state point (p,T) = (1,1) for e = 1.25 and 
a = 0.947, which is a typical liquid state point of the LJ 
system. This was done by proceeding as follows [111]. We 
have previously @, E3] derived the identity 
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(3) 



in which W is the virial, U the potential energy, and 
the angular brackets denote NVT ensemble averages. 
The physical interpretation of 7 (which varies with state 
point) is that n = 3j is the exponent of the inverse- 
power-law potential approximately reproducing the be- 
havior at the state point in question. Combining Eq. ([2]) 
with Eq. §3$) and the simulation results for the fluctua- 
tions of W and U leads to a = 7/2 — 1 = 1.85. Note that 
the p 2 term of Eq. @ is thus negative, which reflects the 
attractive r~ 6 term of the LJ potential. 
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An isomorph is a set of state points with invari- 
ant structure and dynamics in reduced units Q- Via 
appropriate rescaling, however, an isomorph can be 
interpreted differently: as defining a set of different 
LJ pair potentials that give invariant properties at 
the same state point. These are simply two differ- 
ent ways of looking at an invariant Boltzmann factor: 
Using the unit system where ks = f, Eq. (JlJ im- 



,p-^ N )/[f( S )h(p)]) = 



where 



is 



plies that exp 

CX P (-[l//(s)]£i<J«LJ 

the distance between particles i and j. Along an iso- 
morph /(s) is a constant; if we consider the isomorph 
which includes the state point p = T = 1, then given 
the normalization of Eq. ([2]) we have in fact f(s) = 1. 
The shift in interpretation now comes by noticing that 
the same Boltzmann factor is obtained by considering 
a configuration at unity density and unity temperature 
and a family of isomorphic pair potentials v^j(r) = 
Wlj (d _1 / 3 r) /h(d), where we have dropped the tilde from 
positions and replaced p with d to emphasize the shift in 
perspective. These pair potentials are still LJ potentials, 
but with different energy and length parameters. These 
are the potentials plotted in Fig. [TJ . 
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FIG. 3: (a) The AA particle radial distribution function of 
the Kob- Andersen binary Lennard- Jones mixture for a fam- 
ily of isomorphic pair potentials similar to those of Fig. [TJ 

(b) The A particle incoherent intermediate scattering func- 
tion for the same family of potentials as a function of time 
at the wavevector defined from the maximum of g(r). For 
comparison the dotted line shows the WCA prediction [l 21 ] - 

(c) The function x^i^) f° r the A particles for the same pair 
potentials, including the WCA prediction (dotted line). 



The single-component LJ system does not have a broad 
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dynamic range because it cannot be deeply supercooled. 
To test the robustness of the predicted invariance of the 
physics for families of isomorphic pair potentials, we sim- 
ulated also the Kob- Andersen binary LJ (KABLJ) mix- 
ture [13j that is easily supercooled into the highly viscous 
state. For this system the constant a — 1.29 was iden- 
tified from simulations of 1,000 particles at the reference 
state point (p,T) = (1.60,2.00), using again Eq. 
From the function h(p) a family of isomorphic equiva- 
lent pair potentials was generated much like those of Fig. 
[TJ in particular some of them have a vanishingly small 
attraction. 
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FIG. 4: (a) Distribution of forces for the different single- 
component L J potentials of Fig. Q] at the state point (p, T) = 
(1, 1), based on seven independent simulations of the seven 
potentials, (b) Snapshot of the :r-component of the force F x 
on one particle as a function of time. The system simulated is 
defined by e = 1.25 and a = 0.947; F x was subsequently eval- 
uated for the same series of configurations for the six other 
potentials. These figures show that, even though the pair po- 
tentials are quite different, the forces are virtually identical 
(except at the extrema). 



Figure (3][a) shows the AA particle radial distribution 
functions for these different pair potentials. Figure Glib) 
shows the A particle incoherent intermediate scattering 
function; we here included simulations of the Weeks- 
Chandler- Andersen (WCA) system for e = 1.25 and 
(7 = 0.947 (dotted line). Even though the WCA ap- 
proximation has the correct repulsive forces, its physics 
differs considerably from the isomorphic pair potentials 
[l2| . some of which have also (virtually) no attractions. 
We also calculated X±{t)-> a measure of dynamic hetero- 
geneities. The results shown in Fig. [3][c) are more noisy, 
but confirm the predicted invariance of the dynamics for 
the different pair potentials. 



It would require extraordinary abilities to know from 
inspection of Fig. [TJ that these pair potentials have vir- 
tually the same structure and dynamics. The potentials 
have neither the repulsive nor the attractive terms in 
common, so why is it that they have such similar behav- 
ior? The simple answer is that they result in virtually 
the same forces (Fig. 2J. The point is that the force 
on a given particle is a sum of contributions from (pri- 
marily) its nearest neighbors, so plotting merely the pair 
potential can be misleading. Thus, by reference to the 
pair potential alone, one cannot identify separate roles for 
the repulsive and the attractive forces in a many-particle 
system; there simply are no "repulsive" and "attractive" 
forces as such. 

Separate from illuminating in which sense specific roles 
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are played by the repulsive and attractive parts of the 
pair potential 
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the data presented in this paper sup- 
port an interesting approximate scalingprinciple pro- 
posed in 2003 by Young and Andersen [lg. According 
to this, two liquids at the same temperature and den- 
sity with the same particle mass have the same dynamic 
correlation functions if the static pair correlation func- 
tions are the same. Another way of saying it this that 
for fixed density, particle mass, and temperature, the dy- 
namic correlation functions are functionals of the radial 
distribution functions. 
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(P,T) = (1,2) 
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(b) 



(p,T) = (l,4) 
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The above reported simulations focussed for each sys- 
tem on one particular state point. If the potentials in 
Fig. Q] are really to be regarded as equivalent with re- 
spect to structure and dynamics, however, one should 
test also other state points. We have done this briefly 
for the single-component LJ liquid of Fig. [TJ The results 
are shown in Fig. where temperature was increased 
significantly. Clearly, the degree of similarity observed 
at the state point (p,T) — (1,1) is maintained also for 
the other state points (for comparison Fig. 03c) repro- 
duces the {p,T) — (1,1) results from the left panel of 
Fig. Ha)). 



FIG. 5: Radial distribution functions for the potentials of 
Fig. [1] at other state points than the state point (p, T) — 
(1, 1) studied above. For reference we give in each subfigure 
the value of 7 defined in Eq. ©. (a) (p,T) = (1,2); (b) 
(p,T) = (1,4); (c) (p,T) = (1,8) - the 7's reported here are 
those of the state point {p,T) = (1, 1). 
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